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Abstract 

^ Muon and tau neutrinos (u x ) interact with protoneutron star matter only 

via neutral currents and exchange energy with the stellar gas predominantly by 
neutrino-electron scattering and neutrino-pair processes. In contrast, electron 
neutrinos and antineutrinos (v e and v e ) are frequently absorbed and produced 
in charged-current mediated reactions with nucleons. Therefore the emergent v e 

Q and v e originate from layers with lower temperatures further out in the star and 

CO are emitted with much lower characteristic spectral temperatures. In addition, 

a major contribution to the v e and v e opacities is due to absorptions, while the 

Q\ opacity of v x is strongly dominated by scattering reactions with nucleons and 

nuclei in which the v x energy is (essentially) conserved. Therefore the v x dis- 
Qh tribution is nearly isotropic when v x decouple energetically and their outward 

q diffusion is slowed down. In a generalized form to include this effect, the Stefan- 

j^j Boltzmann Radiation Law can account for both the facts that v e (v e ) and v x 

^ emerge from the star with similar luminosities but with very different character- 

istic spectral temperatures. Simple analytical expressions to estimate the effect 
are given. If, as recently argued, even at densities significantly below nuclear 
matter density neutral-current scatterings were associated with considerable en- 
ergy transfer between neutrino and target particle, one might expect spectral 
temperatures of v x much closer to those of v e and v e . This is of relevance for the 
detection of neutrino signals from supernovae. 



*On leave from Max-Planck- Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, D-85740 Garching, 
Germany. 
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1 Introduction 



The transport of v e and v e on the one hand and of heavy-lepton neutrinos v^, z^, u T , 
and v T on the other is different in protoneutron star matter. Heavy-lepton neutrinos 
can be handled as one species v x , because for vanishing or small (of the order of 
few eV) neutrino masses differences of their behavior are insignificant. Muons and 
tau mesons are not present in the outer regions of protoneutron stars; the maximum 
temperatures T in newly formed neutron stars are several 10 MeV only and at densities 
P & Pnuc (Pnuc ~ 2.5 • 10 14 g/cm 3 is the nuclear matter density) the electron chemical 
potential is usually below 100 MeV, too low to allow for a significant production of 
muons. Therefore electron-type neutrinos interact via neutral and charged currents 
with the particles of the stellar gas, while heavy-lepton neutrinos can react only through 
neutral currents. This difference affects neutrino-electron scattering and neutrino-pair 
processes, which are the most important reactions by which v x exchange energy with 
the stellar gas. The main effect on the transport, however, results from absorptions 
of v e and v e onto n and p, respectively. These reactions (and the inverse neutrino- 
producing processes) enforce reactive and thermal equilibrium between v e and v e and 
the stellar gas until further out in the star, to lower densities and temperatures, where 
the weaker 7vt>-pair production processes are already frozen out for all neutrinos and 
where v x are therefore energetically decoupled from the background. 

u e , V> e , and v x with the same energy scatter off n or p with the same cross sec- 
tions. Since neutrino-electron scattering and neutrino-pair annihilation yield only a 
minor contribution (< 10%) to the total transport opacity, the latter is dominated 
by neutrino- nucleon scatterings in case of v x . In contrast, neutrino absorptions onto 
baryons provide a sizable fraction [k, 30-80%) of the transport opacity in case of v e 
and i> e . This leads to different characteristics of the particle transport through the 
semi-transparent outer stellar layers, where the spectral distributions of the emitted 
neutrinos are determined. 

v e and v e emerge from regions with lower temperatures further out in the star and 
are emitted with much smaller characteristic spectral temperatures than v x . Never- 
theless, detailed models and transport simulations [1] show that all neutrino types are 
emitted with very similar luminosities. Unless v e and i> e emission from accretion onto 
the protoneutron star plays an important role the luminosities are equal to much bet- 
ter than a factor of 2. (For a review of the model results and more details, see [2].) 
These findings will be discussed and explained by the different transport properties of 
v e and v x in the presented paper. Simple analytical relations will be presented that 
allow one to estimate the size of the considered effect (Sect. 6). The implications and 
consequences will be addressed in Sect. 7. 
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2 The neutrino "spheres" 



The differences of the interactions as described above lead to different relative loca- 
tions of the "energy spheres" and "transport spheres" for neutrinos of different types. 
The "energy sphere" is defined as the layer in the star external to which a neutrino 
of representative energy, before leaving the star, interacts only once by a reaction 
with energy transfer to or from the target particle (optical depth for reactions with 
energy exchange ~ 1). The transport sphere is the corresponding layer with optical 
depth about unity for all reactions, including those which leave the neutrino energy 
(essentially) unchanged ("conservative processes"). 

Since a significant fraction of the total transport opacity of v e and v e is due to 
absorptions onto nucleons, this process, however, is absent in case of v x , the energy 
and transport spheres of v e and v e lie very close together and further out at lower 
temperature and density. For v x both spheres are spatially clearly separated by the 
strong dominance of conservative scattering reactions in the transport opacity. Typi- 
cally, from inside out the v x energy sphere, the v x scattering sphere, the v e (i> e ) energy 
sphere, and the v e (p e ) scattering sphere follow in the direction of decreasing temper- 
ature and density. Dependent on the ratio of the mean spectral energies of v e and v e 
relative to v x and caused by the variation of absorption and scattering cross sections 
with roughly the square of the neutrino energy, the radial locations of v x scattering 
sphere and v e (i> e ) energy sphere may also follow in inverse order. For steep gradients 
of T and p at the protoneutron star surface, the differences of the radii may be small, 
while the density p and temperature T change appreciably. Initially, the transport 
character of v e is very similar to that of v e . As the deleptonization goes on and the 
protoneutron star matter gets increasingly neutron-rich, however, it becomes more and 
more similar to that of v x . 

3 The "paradox" 

As the protoneutron star deleptonizes and the electron chemical potential p, e and the 
electron fraction Y e decrease, the chemical potential of electron neutrinos approaches 
zero: p, Ue = p e + p p — (t n — > 0. By this time the total energy fluxes (luminosities) 
of all neutrinos (v e , v e) and v x ) become very similar in all detailed numerical models, 
L„ e L 9e rs L„ x . Nevertheless, v e (and P e ) leave the star with significantly lower 
average energies than v x , typically (e^) r; (2 ... 2.5) • (e„ e ) r; (1.3 ... 1.7) • (e^). 

To state the paradox in the above situation more clearly, let us assume that neutri- 
nos stay roughly in thermodynamical equilibrium with the stellar gas until they decou- 
ple energetically at their corresponding energy spheres (narrow decoupling layers) at 
radii R eVi (z^ = v e , v e) v x ). This means that to zeroth order of accuracy the neutrino 
spectra can be represented by Fermi-Dirac distributions according to local equilibrium, 
the neutrino temperature being equal to the local gas temperature, T Vi = T(R EjVi ). 
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Moreover, since at "late" times the chemical potential of electron neutrinos is small 
in a large part of the protoneutron star, /i Ue = —/jp e ~ 0, and \x Vx = is true for the 
v x chemical potential at all times, the emission of neutrinos from their corresponding 
energy spheres can roughly be considered as blackbody-like with (e Vi ) oc T(R £:Vi ). v x 
decouple at a somewhat smaller radius, but higher temperature, which is reflected by 
T(R e ^ x )/T(R e ^J <: 2. From the Stefan-Boltzmann Radiation Law one calculates the 
total neutrino luminosity according to L v . oc R^ v .T 4 (R e ^.). Since differences of the 
energy-sphere radii of different neutrinos are very small compared to the neutron star 
radius (because of the steep density and temperature gradients close to the surface), 
i.e. R t , Ve ~ Rt,v e ~ Re,v x ~ Rns, one ends up with the conclusion that L„ x ^ 16 • L„ e . 
This is in strong contradiction to the results of self-consistent model calculations ac- 
cording to which the luminosities are equal within a factor of less than about 2. How 
can this be understood? 

4 Protoneutron star cooling 

The emergent neutrino flux is certainly not a perfect equilibrium emission from a well 
defined neutrino sphere but is much more complex with a non-equilibrium spectral 
shape and with neutrinos of different energies decoupling from the stellar gas at dif- 
ferent radii and at different thermodynamical conditions, simply because the neutrino 
opacity varies roughly with the square of the neutrino energy. This might serve as 
an "explanation" of the described puzzle, but does not actually yield insight into the 
direct reason why the neutrino luminosities can be very similar despite of the spectral 
dissimilarities. In fact, there is a simple way to understand qualitatively the physical 
mechanism which is responsible for the observed result. 

Firstly, note that in the surface layers of the protoneutron star, even somewhat 
inside of the neutrino energy spheres, the total neutrino energy flux (as the sum of the 
luminosities of all neutrinos and with general relativistic redshift being disregarded) 
is conserved and given by the energy-loss rate of the dense core region of the star. 
External to the corresponding neutrino-energy sphere the total energy flux of each 
individual neutrino type is constant. The heat capacity of the dense core being very 
large, it is the energy loss from the central part of the star that determines the Kelvin- 
Helmholtz time scale of neutrino cooling. During the Kelvin-Helmholtz phase the 
diffusion of neutrinos out of the very opaque regions governs the evolution. Due to the 
long diffusion time scales (order of seconds) the cooling proceeds quasi-stationarily with 
very slow changes of thermodynamical quantities and neutrino emission parameters. 
In contrast, the surface layers heat up or cool down quickly when energy is deposited 
or lost via neutrinos. Within a time scale of only several 10 ms the temperature and 
density structure re-adjusts by expansion or contraction to assure that the neutrino- 
energy flow diffusing out from deeper regions is just transported through the surface 
layers. 
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Note, secondly, that for analogous reasons there is a conservation of the total net 
lepton-number fluxes (total number flux of neutrinos minus antineutrinos of each flavor) 
in the surface layers, but no conservation of the total number fluxes of individual 
neutrinos. Under the constraints that the total net lepton-number flux is conserved 
and the Pauli exclusion principle is obeyed, neutrinos are absorbed and re-emitted 
and redistributed in phase space. This is completely analogous to photons when they 
diffuse to the surface of a star with conserved total energy flux but with a down- 
grading in energy space as the temperature decreases towards the surface. The constant 
total energy flux is transported by an increasing number of neutrinos or photons. 
The redistribution in energy space determines the spectrum with which neutrinos are 
finally emitted from the neutrino-energy sphere, external to which energy changing 
interactions become very rare. 



5 Fundamentals of radiative transfer 

Considering neutrino diffusion through the background medium the neutrino energy 
flux is given by 



D- V 



4 V 4 fs(0) 



(he) 



-~4~ -3^3(0) T 3 VT, (1) 
3 (he) 



where D is the diffusion constant, D = A t c/3, with the total transport mean free 
path A t (as determined from all neutrino reactions). T is the gas temperature. Since 
neutrinos are assumed to be in thermal equilibrium with the stellar gas the considered 
transport approximation is called "equilibrium diffusion". ^3(0) denotes the third 
Fermi integral for relativistic particles with degeneracy parameter 77 = [ijT = 0, the 
general definition being T n (rj) = dx x n [1 + exp (x — r/)] -1 . One can express the 
energy flux in terms of the "flux factor" f and the neutrino energy density E by 

Air 

F = fcE = ie— -3^3(0) , (2) 
(he) 

when f is defined as the normalized flux, 

f = / dnni(n)/ / dfti(n) . (3) 

In spherical symmetry / as the first angular moment of the radiation intensity 1(H) 
gives the mean angle cosine of the radiation relative to the radial direction r/r = n: 

f = / . n = (cosfl) n . (4) 

Combining Eq. (1) and Eq. (2) we find 

4 VT 

f = --.V-. (5) 
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With Eq. (2) the neutrino luminosity of the protoneutron star can be expressed in 
terms of the conditions at R e , Vi as 

L Vi = 4nRl Vi -f Ui (R e ^-c — -T 4 (R e , Vi )-F 3 (0) 

oc Rl Vi .f Vi {R e , Vi )-T\R e , Vi ) . (6) 

If one assumes diffusion to be a suitable description until the energy flux decouples at 
the neutrino-energy sphere, f Vi (R e ^ Vi ) is given by Eq. (5): 



fvi{Re, Vi ) — — ~ • \, Vi {Re, Vi ) 



3 \ e '"« / ^ r 



(7) 



Re, Vi 

In analogy to the photon case one can define an effective temperature T* s by 

f Vi = k = a*(r;f) 4 (8) 

with the radiation constant a* for neutrinos given by 

a* is the neutrino-density constant corresponding to the radiation-density constant 
a. The effective temperature T£ ff is related to the characteristic spectral temperature 
T(R^ V .) by comparison of Eq. (8) and Eq. (9) with Eq. (6): 

Ttf = (4U,(R^)) 1/4 -T(R^) . (10) 

For causality reasons (i.e. F < Ec) and by Eq. (4) / is physically limited to values 
/ < 1. Close to and inside of the neutrino transport sphere one typically finds / ^ 1/4 
[3, 4]. Therefore Eq. (10) implies that T£ ff ^ T(R e ^.). Note that it is the characteristic 
spectral temperature T u . = T(R e ^ u .) that is connected with the average energy (e Ui ) of 
the emitted neutrinos according to 

(e ) = T -^21 (11) 
if one assumes a vanishing spectral degeneracy parameter, r\ Vi = 0. 



6 The explanation 

Eq. (6) and Eq. (7) allow us now to resolve the paradox of Sect. 3. From Eq. (6) one 
learns that the luminosity is not just proportional to the product R^.T 4 (Re,^) , but 
in a general treatment the flux factor f(R^ v .) appears. This is usually disregarded 
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and f(R €;Vi ) is replaced by the value one obtains for a (black) spherical surface that 
radiates isotropically into vacuum: 

/ = P""//- 1 , 11 " = 3 (12) 

("vacuum approximation") 1 . The implicit assumption f Vx (R eyVx ) = f Ve (R etVe ) is not 
justified in the case of v x and v e due to the differences of their transport behavior as 
discussed in Sects. 1 and 2. In detailed neutrino transport calculations with Monte 
Carlo methods [3, 4] values of f Vi {R\,v?) ~ 1/4 are obtained for all types of neutrinos 
at their corresponding transport spheres but usually not at their energy spheres. 

At the neutrino-energy spheres (index e) and transport spheres (index t) defined 
by optical depths of order unity, r t)e ~ 1, the approximate relations hold: 

T~t,e ~ ^t,e ' Rt,e ~ 1 M,e ~ Rt,e ■ (13) 

A t and A e are the transport mean free path and the mean free path for reactions with 
energy exchange, respectively, and K t ,e = A^ 1 are the corresponding opacities 2 . The 
v x transport mean free path \, Vx is much smaller than the mean free path for energy 
exchange, because the opacity is dominated by conservative scattering reactions. Thus, 
< A e ,„ x (r), therefore \ Vx {R t , Vx ) « R t ,„ x « R ns (Eq. (13)) and \ e ^{R e ^) « 
Re,u x ~ Rns, but X t , Vx {Re,u x ) < Rns- In contrast, for v e (or v e ) one has \ Ve {r) ~ K,u e (r) 
and therefore A t:I/e (R e ,u e ) ~ Rns- Putting all together one infers from Eq. (7) 

fMr.) « U e (Re,u e ) ■ (14) 

Using this in Eq. (6) finally explains, why L Vx L Ve can be true although T„ x 3> 
T Vb . This is simply a consequence of the isotropization of the v x distribution due to 
frequent (iso-energetic) scatterings off nucleons, which slows down the v x diffusion and 
suppresses the v x flux from the v x energy sphere. Although v x decouple energetically 
deeper inside the protoneutron star and have a more energetic spectrum that reflects 
the higher temperature at the v x energy sphere, this leads to a v x luminosity similar 
to that of v e . 

To obtain a numerical estimate of the corresponding effect represented by Eq. (14), 
Eq. (7) can be evaluated for typical v x and v e energies with the relevant processes 
included. With the typical number (e Ux ) / (e Ue ) ~ 2 ... 2.5 (Sect. 3) and with taking into 



■^ote that this is actually not a consistent definition, because one takes into account in the numer- 
ator only the part of the radiation that is directed outwards. For isotropic radiation the integration 
from —1 to 1 would yield / = 0. Moreover, moving infinitesimally outward from the radiating surface 
one obtains / = \. Thus / is discontinuous. For an introductory description in the context of photon 
transport, see, e.g., [5]. 

2 Being precise, the determination of the neutrino-energy sphere requires the use of the effective 
mean free path for energy exchange, A| ff = \/\ A c , instead of A c . A| ff represents the average displace- 
ment between two reactions with energy exchange (see, e.g., [5, 6, 3]). The presented arguments are 
qualitatively not affected by this. 
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account that the scattering opacity yields a fraction of 20%-70% of the total transport 
opacity of v e (or v e ) (Sect. 1) and that the spectral averages of the scattering and 
absorption cross sections vary roughly with (e 2 ), one estimates for the relation of the 
cross sections: a t ^ x (R etVx ) m (1 ... A) ■ a t ^ e {R ttVi! ) . The cross sections become comparable 
when absorptions dominate scattering in the v e opacity. From this follows 



Since R e , Vx < Re, Ve one always has p(R e , Ve ) / p(R e ^ x ) < 1 for the ratio of the densi- 
ties at the v e and v x energy spheres. For steep density profiles the density ratio can 
become quite small, which is the case at late times during the protoneutron star cool- 
ing. At early times the neutronization of the star is not so strong and v e absorptions 
onto n less frequent. Therefore the total transport cross section of v e and the ratio 
Vt,v e ( R e,v e )/vt,vA R e,vJ are smaller. 

7 Discussion and implications 

The result of detailed protoneutron star cooling models that the luminosities of emer- 
gent v x and v e (p e ) are similar although the spectra are very different, can be un- 
derstood as a consequence of the differences of the interaction of v x and v e (p e ) with 
protoneutron star matter. v x interact by neutral currents only, while v e and v e react 
via neutral and charged currents. The v x transport opacity is therefore dominated by 
(essentially) iso-energetic neutrino-n and neutrino-p scatterings, whereas absorptions 
onto nucleons yield a sizable contribution to the opacity of v e and i> e . Due to the 
stronger coupling to the stellar matter v e and v e emerge from the star further out, at 
lower temperatures and densities than v x . All their interactions with the stellar gas 
cease at around the same radius. In contrast, when v x decouple energetically from the 
stellar gas at their energy sphere they still experience frequent scatterings off nucleons 
(and nuclei) on their way out. This "back-scattering" causes an isotropization of the 
v x distribution even exterior to the v x energy sphere and slows down the diffusion of 
v x to the stellar surface. The corresponding suppression of the v x luminosity leads to a 
significant discrepancy between the effective temperature (oc L]/ 4 ) and the charac- 
teristic spectral temperature T Vx (oc (e„ x )). The emergent v x spectrum remains nearly 
unchanged outside of the v x energy sphere and T Vx therefore reflects the temperature 
in a deeper layer of the star than T Ve or T Pe . 

Instead of absorbing the effects of scattering-dominated v x transport into the effec- 
tive temperature one can explicitly account for the different angular distributions of v e 
(i> e ) and v x by using a generalized form of the Stefan-Boltzmann Radiation Law that 
includes the flux factors f Vi (R etVi ). Three parameters, f Vi (R^ m ), T Vi = T(R e ^ i ), and 





P( R e,u x ) C"t,^(-Re^J 




(15) 
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R^ Vi fa R ns characterize the physical dependences of the luminosities. Considering the 
transport mean free paths of v e (p e ) and v x at their corresponding energy spheres a 
simple analytical way to estimate the influence of the discussed effects on the properties 
of the emergent neutrino fluxes is given in Eq. (13)-Eq. (15). 

A potential difference between the characteristic spectral temperature and the ef- 
fective temperature according to Eq. (10) has been neglected in all evaluations of the 
SN 1987A neutrino signal that were based on simple cooling models of the nascent neu- 
tron star (e.g., [7]). In case of v e this difference is quite pronounced at late times when 
the protoneutron star matter has become very n-rich and the v e and v x transport are 
more similar than during the early phase of the cooling evolution. In time-dependent 
analyses of the neutrino data the effect can therefore not be ignored. 

Deviations from an isotropic radiation field at the emitting surface account for 
the phenomenon of "limb darkening" which is strong when / tends to values close to 
1, indicating an outward peaked angular distribution of the radiation. In case of an 
extended, radiating surface the limb darkening effect in combination with the char- 
acteristic radiation temperature in principle allows one to determine the temperature 
structure as a function of depth below the surface (Eq. (7)). A reliable value of the 
factor f(R £tV .) is also crucial to recover neutron star radii from measured neutrino 
luminosities and characteristic spectral temperatures (Eq. (6); see also [2, 3]). 

The addressed differences of v x and v e (u e ) transport and the associated conse- 
quences for neutrino fluxes and spectra were based on the assumption that the en- 
ergy exchange between neutrinos and target particles in neutrino-nucleon scatterings 
is negligible. However, as recently argued [8], this need not be true because neutri- 
nos interact with nucleons in a dense plasma and not with isolated targets in a very 
dilute medium. For typical protoneutron star temperatures even at densities as low 
as 10 13 g/cm 3 nucleon-nucleon collisions could allow for significant energy transfer be- 
tween nucleons and neutrinos that is not possible in interactions with isolated nucleons 
due to the smallness of nucleon recoils. In the context of the presented paper this has 
the important implication that, if different neutrinos are not affected in the same way 
or not equally strongly, the differences between the spectra of v x and v e or v e could 
be significantly changed. In the considered "standard" case of conservative neutrino- 
nucleon scatterings the neutrino-energy spheres, in particular the energy sphere of v x , 
move inward to very high densities as the temperatures in the protoneutron star de- 
crease along with the cooling process and the decrease of the neutrino opacities due to 
Pauli blocking in the Fermion phase spaces. Including the effects of nucleon-nucleon 
interactions in the dense medium might lead to more similar v x and v e spectra when 
the discrepancy between the v x transport mean free path and the mean free path for 
energy exchange shrinks. A quantitative analysis of this possibility needs more theo- 
retical work on the cross sections; incorporation into the numerical modelling of the 
neutrino cooling of protoneutron stars might yield very interesting results. It is well 
possible that detections of neutrinos of all flavors from a future supernova, e.g. in 
the Superkamiokande and SNO experiments, might allow to determine the size of the 
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discussed effect by a comparison of the v e and v x spectra. 
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